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I. INTRODUCTION 



The use of light and lasers in medicine have increased manifold in the past decade. Enormous efforts have been 
devoted to the development of new diagnostic techniques such as NIR imaging 0] and flourescence spectroscopy of 
tissues 0], as well as therapeutic uses such as photothermal coagulation and Photo-dynamic therapy ||. Tissues 
are highly scattering media and an accurate knowledge of light transport parameters of the tissue is indispensible 
to describe the propagation of light in these media. Several techniques have been developed to measure the light 
transport parameters in tissues B. Apart from these, the use of the Coherent Backscattered peak from tissues to 
estimate the transport parameters has also been suggested 

The phenomenon of Coherent Backscattering (CBS) of light by random media has attracted considerable attention 
since 1985 |]||, when connection was first made between elastic multiple scattering and Anderson localization. This 
phenomenon also termed weak localization, shows up as a sharp peak in the backscattered direction within a narrow 
cone of angles. The angular width of the CBS peak is mainly determined by the transport length of light in the 
medium and the shape of the peak is slightly modified by the absorption. This technique was used to characterize 
the transport lengths of samples used in experiments on strong localization of light ]1C[ ] . Etemad et. at, jllj were the 
first to study the effects of absorption in random media on the CBS peak. 

There have been some attempts to study the CBS peak from biological tissues and obtain the transport properties 
|5|-@] . These studies, however, have been inaccurate because of a misinterpretation in applying Akkerman's ex- 
pression ]l2] | derived for conservative, semi-infinite media to absorptive media. Hence a systematic study of the CBS 
peaks from biological tissues has not yet been made. Since this is a non-contact method which can be used for online 
in-vivo measurements also, it is important that a systematic evaluation of this technique for the measurement of 
transport parameters be carried out. Eddowes et. ai, |]l4| have suggested that accurate absorption information may 
be obtained by using Monte-Carlo simulations to fit experimental measurements and have developed Monte-Carlo 
based routines to identify the optical coefficients for a given CBS angular profile. However the measurement of the 
absorption parameter is dependent on both the transport length and the limitations of the experimental setup, as 
will be explained, and accurate absorption information can be obtained in only limited ranges of the absorption and 
transport coefficients. In reports that have appeared so far, researchers have concentrated on obtaining the transport 
and inelastic lengths of tissues from the shape of the CBS angular profile. The possibility of using the intensity 
information for estimating the anisotropy factor^) of the scatterers has not been examined. 

In this paper, we present the results of our investigations on the suitability of using CBS for measuring transport 
parameters of biological media. Milk of different concentrations have been used as tissue phantom. It is known that 
the light transport properties of milk closely resemble those of tissues and milk has been used to model light transport 
in brain tissue p"5| , to study optical imaging Jl6[ , photon-density waves , and for investigating the propagation of 
short laser pulses through a scattering medium [[18|. To investigate the effects of absorption on CBS, an absorbing 
dye with known extinction coefficient was added to milk and the CBS from such media were studied. We also exam- 
ine the possibility of determining the anisotropy factor (g) from the the single scattering contribution to the diffuse 
background. 



II. COHERENT BACKSCATTERING OF LIGHT 



Coherent Backscattering (CBS) of light occurs in all disordered media and is the only major surviving interference 
effect. When a beam of light is incident on a random medium, there exist partial waves traversing every possible path 
in the medium. The CBS effect arises from the constructive interference of any partial wave with its time reversed 
counterpart in the medium. In exactly the backscattered direction, both these two waves have the same phase and 
constructive interference results. Away from the backscattered direction, the counterpropagating paths develop a 
phase difference depending on the relative positions of the first and last scattering events in the medium. For the 
ensemble of all possible light paths, these phases will randomize and the reflection is enhanced within a narrow cone 
in the backward direction with an angular width of the order of X/lt where A is the wavelength of light and It is the 
transport length in the medium. This peak shows up only after the ensemble averaging over the large scale sample 
specific fluctuations (speckle) that originate from the random medium fl9|| . 
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The CBS intensity can be described in terms of three contributions. The total normalized angular intensity is 
described |^(| by 

r<n n lc(e,L) +ll (e) + ls {6) 
7/(0) + 7s (0) 

where j c , ji and "f s are the bistatic coefficients of the coherent, diffuse and single scattering contributions respectively. 
9 is the angle of scattering and L represents the set of transport lengths in the medium. The set of relevant transport 
parameters are the mean scattering length(Z s = l/(pa s )) defined as the reciprocal of the average of the product of 
the single particle scattering cross-section (cr s ) and the density of scatterers (p)(the probability of a photon remain- 
ing unscattered after traversing a distance z in the medium is exp(— z/l s )), the inelastic length (/,•) defined as the 
reciprocal of the absorption coefficient (a) (transmitted intensity It(z) — Iq exp(— z/li) in the absence of scattering), 
the anisotropy factor (g = (cos 9)) which is defined as the average of the cosine of the scattering angle for a single 
scattering event and the mean transport length (l t ) which is a measure of the distance in which, the direction of 
the photon's motion becomes uncorrelated with its initial direction and is related as It — Z s /(1 — g). The bistatic 
coefficients for the diffuse and coherent intensities are traditionally described by summing up the ladder diagrams and 
the most crossed diagrams respectively in a systematic perturbation of the intensity propagator The bistatic 

coefficient for the single scattering contribution for isotropic scatterers assuming normal incidence is given by p2| 
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where A is the area of the target, no is the density of scatterers, a s is the single particle total scattering cross-section, 
o a is the absorption cross section, a is the albedo defined as the ratio of the total scattering cross-section to the 
extinction cross-section[er s /(er s + cr a )] and b{= no(a s + cr a )L) is the optical thickness of the slab. However for the case 
of anisotropic scatterers, because it is the backscattering coefficient which would contribute to the measurement, a 
would be modified as 

a = (3) 

<Js+<?a 

where ab is the single particle back-scattering coefficient. 

When the single scattering contribution is supressed by using an isolator Jl9"| and the helicity preserving channel 
for circularly polarized light is detected from time invariant random media, the enhancement is exactly 2 in the 
weak disorder regime pp| ]. However if the detection is carried out in the linear polarization preserving channel, the 
enhancement factor will always be less than 2. This is due to the presence of single scattered events which do not 
contribute to the enhancement in the back-scattering. In earlier work J^-Q], the contribution of single scattering was 
ignored. Though Eddowes et. al. [[u} suggest that the experiment must be performed with circularly polarized light, 
it has not been examined whether the single scattered contribution can be used to estimate the anisotropy factor. 
The dependence of a (normalized to unity for isotropic lossless scatterers) with the anisotropy factor g (obtained by 
varying the scatterer size), is shown in figure-1. The cross-sections and the anisotropy factor were computed using Mie 
scattering theory for spherical particles and a program based on the BHMIE code given in |22j ]. A value of 1.33 for 
the refractive index of the medium (water) and 1.41 for the refractive index of the particles (typical of tissues Q|) is 
used. Note the presence of resonances for monodisperse particles in figure-la. In figure-lb, the cross-sections and the 
g factors were then averaged over a particle size dispersion of 30% assuming a Gaussian distribution in order to mimic 
the experimental situation where the particles are polydisperse . Slightly different values of the refractive indices 
did not lead to much change, as the averaging process washed all resonances present. As one can notice, the single 
scattering contribution decreases sharply with the increase in anisotropy and becomes very small for large anisotropy 
factors. For g>0.6, the single scattering contribution from lossless scatterers is barely 2% of the total back-scattered 
intensity and would hardly be measurable. This unfortunately implies that the anisotropy factor of biological tissues 
(typically in the range 0.8-0.95 H|) would hardly be measurable by this technique. 



The CBS peak shape is reasonably well described in the diffusion approximation |21| , |12fl . In this work, the expres- 
sions derived in the diffusion approximation for an incident plane wave on a slab of isotropic scatterers and with a 
finite thickness, have been used [^l| (See Appendix for the expressions). By replacing the mean scattering length (l s ) 
in these expressions by the mean transport length (Z t ), these expressions have been shown to be valid for anisotropic 
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scatterers as well ]24| . 

It is essential to observe the effects of absorption on the CBS peak. The main effect of absorption is a rounding of 
the central cusp of the peak in an angular range 9 a ~ \j \Jl t li. This corresponds to the extinction of the longer path 
lengths greater than y/lth due to the absorption. The other effect of absorption is to trivially reduce the enhancement 
factor because the single scattering events become relatively more important as absorption increases. The CBS peak 
is also rounded by the convolution with the instrumental response function and unless the angular range 9 a of the 
rounding due to absorption is larger than the instrumental resolution, the absorption would not be measurable accu- 
rately. For the case of 6 a = 0.l8mrad our experimental resolution, A = 0.6328/im and It = 670/im a minimum value 
absorption coefficient a ~ 0.6cm -1 can just be discerned. Though much smaller absorption coefficients in principle 
can be measured by improving the angular resolution, even assuming ideal optics, the resolution eventually is limited 
by diffraction due to the finite beam size. For our beam size of 5.5 mm, a beam divergence (64 = X/nR where R is 
the beam radius) of about 75fj,rad is present, yielding a minimum measurable a of about 0.1cm -1 at a It of 700/^m. 
Trying to increase the beam size further is impractical due to problems of inhomogeneity of tissue samples. Hence 
accurate information would be possible only for samples with moderate absorption (a > 1cm -1 ). However this is by 
no means a fundamental limit of the technique and depends on the transport length of the sample and the resolution 
of the experimental setup. 



III. EXPERIMENTAL SETUP AND PROCEDURES 

The experimental setup used is schematically depicted in figure-2. A 5mW He-Ne laser at 632. 8nm was expanded to 
a beam diameter of 5.5 mm and collimated to diffraction limit. The collimation was checked by shear interferometry. 
Following the standard practice, the CBS light was viewed through a non-polarizing 50-50 beam-splitter with a small 
wedge. A computer controlled CCD array (752 X 244 pixels, model EDC-1000HR, Elcctrim corporation, U.S. A) was 
placed at the focal plane of a positive lens with a focal length of 20 cm, to analyze the angular peak. Each pixel on 
the CCD in this configuration corresponds to 0.06mrad. The laser beam was linearly polarized by placing a polarizer 
just before the beam-splitter. Another polarizer was placed in front of the CCD to record images in the polarization 
preserving channel. An aperture slightly larger than the input beam was placed just behind the beam-splitter to avoid 
ghost images from the AR coated surface of the beam-splitter. The free beam going through the beam-splitter was 
carefully damped using ND filters. 

The setup was aligned by placing a mirror in place of the sample and making the reflected beam go back into the 
laser. The beam transmitted by the beamsplitter and focussed by the lens was scanned by the CCD at the focal plane 
and the intensity profile so obtained was used as the resolution curve characterizing the system response. This was 
well described by a gaussian having a FWHM of 0.17 mrad. The tissue phantom used was commercially available 
boiled and skimmed milk. The concentration of the milk was varied by adding distilled water. To prepare samples 
with different absorptions, the milk was doped with known concentrations of methylene blue dye. Methylene blue 
is a water soluble dye having a broad absorption band in the red region with two peaks at 610nm and 664nm The 
extinction coefficient of the methylene blue solutions which were added to the milk, was measured by a Shimadzu 
spectrophotometer. The dye-doped colloids thus had a well determined inelastic length. The milk did not change 
its scattering properties measurably in a time of about 3 hours. This was confirmed by recording the CBS profiles 
at different times and comparing them. All the experiments were hence carried out on the same milk sample within 
this time frame. The tissue phantom was taken in a cuvette of 10mm path length. It was placed slightly tilted to the 
incident beam so that the specular reflection was well away from the backscattered direction. The Brownian motion 
of the milk particles caused an ensemble average over the speckle and sharp symmetric peaks were observed. CBS 
intensities for samples of different scatterer concentration and absorption were recorded. 

The CBS peaks from milk suspensions were very sharp (FWHM ~ lmrad) and consequently affected by the finite 
instrumental response. Hence in order to fit the parameters, the theoretical profiles were first convolved with the 
system resolution curve using FFT routines. The resulting curve was then least squares fitted to the experimental 
points treating l t ,U and 7,5 as parameters to be fitted using standard NAG library routines. Deconvolution of the 
experimental points was avoided as this was found to be a noisy process. 

To supplement the CBS measurements, independent measurements of the transport parameters were made by 
measuring the transmission of light through milk solution as a function of the optical thickness or the concentration 
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of milk. The milk was taken in a cuvette of 5mm path length. A 2mW He-Ne laser was used as the source. The light 
exiting at the back of the cuvette was imaged onto a photodiode, whose signal was detected by a lockin amplifier 
controlled by a computer. A separate photodiode was used to simultaneously record the laser intensity fluctuations. 
We have also undertaken angle-resolved scattering measurements for estimating the value of the anisotropy factor at 
the same wavelength. The sample was taken in a 1.0 mm cylinderical cuvette and is diluted enough, so as to get into 
the single scattering regime. The transmission of the sample was coupled to the PMT through an optical fibre kept 
on a rotating stage such that the distance from the sample to the fibre tip does not change and is read by a digital 
storage oscilloscope. A reference detector monitored the laser intensity fluctuations. Measurements were avoided in 
the lower angle region upto 8 degrees. 

IV. RESULTS 

Figure-3 shows the CBS profiles obtained from two different concentrations of milk. The solid lines shown are the 
best fits to the corresponding experimental points. As expected the CBS profiles became narrower with increasing 
transport length. The enhancement factor for the sample with longer transport length is considerably reduced due to 
the effect of convolution with the system response. The fits are quite good with a figure of merit (\ 2 /N where N is 
the number of degrees of freedom for the fit) of about 1.2. The variation of the inverse of the transport length with 
the concentration is shown in the inset of figure-3. The l^ 1 increases linearly with concentration as expected. The 
transport length is well defined and the numerical routines used quickly converge to the correct value to an accuracy 
of ±5% in a few iterations. However for the plain milk solutions any absorption coefficient less than 1.5cto _1 could be 
fitted with marginal changes in the transport length. This is because the intensity measurement accuracy is limited 
to about 3% by the CCD noise and the changes in the profile caused by small amounts of absorption are of the same 
level. 

Figure-4 shows the CBS intensity obtained from dye-doped milk with a pre-determined absorption coefficient of 
5 cm" 1 . The fitted curve is for an absorption coefficient of (5.2 ± 0.2)cto _1 and the elastic transport lengths are 
within reasonable variation for the pure milk (l t = 625 ± 30/ito) and the absorbing sample (It — 670 ± 30/ito). At 
higher absorption, the entire peak appears broadened because the angular range of the rounded cusp is 9 a ~ O.Gmrad 
while the peak width itself is of the same order. Also the height of the peak seems to be relatively unaffected by the 
absorption even though the weight of single scattering is enhanced. This counter-intuitive behaviour occurs because 
of the decreased effect of the instrumental response on the broader peak compared to that of the lower absorption 
narrow peak. This is clearly brought out in figure-5 where the theoretical shape (assuming same 7 S as from experi- 
ment) for an infinite resolution is plotted. The high absorption peak clearly shows a lower enhancement. Thus the 
finite resolution and the relative weight of single scattering act in opposite directions on the enhancement factor with 
increasing absorption. The inset in Figure-4 shows the CBS profile and the theoretical fit from a sample with higher 
absorption (a = 10cm -1 ). The backscattered intensity level was small necessitating large exposure times, increasing 
the CCD noise. However the theoretical fit again correctly yielded the absorption coefficient within an accuracy of 
±10%. We studied the CBS in dye doped milk having absorption coefficients upto 30cm _1 . However, the noise 
at these absorption levels was high (around 30 % of the peak intensity) and would appear to limit the use of this 
technique to measure larger absorption coefficients. Anyhow at higher levels of absorption, the expressions for the 
CBS peak shape derived within the diffusion approximation will no longer be valid. 

Using the best fitted values of j s , the enhancement factors of the deconvolved curves if the plain milk- water solu- 
tions (no dye) are found to be in the range of 1.82 to 1.86. This reduction in the enhancement factor is lesser than 
the case of isotropic scatterers in a finite slab of same thickness with the same It and k for which the theoretical 
enhancement factors are between 1.76 to 1.80. These values are in agreement with those reported by Wolf et. al. [ p5| , 
but as observed by them the reduction in the enhancement factor appears much too large for such forward scattering 
media. Other processes such as recurrent multiple scattering pfj ] which reduce the enhancement are not important 
at such large mean free paths present here. The lowering of the enhancement factor could also be caused by the 
finite beam effects and the gaussian intensity profile of the laser beam. However we do not believe this to be the 
reason as the beam size is reasonably large (5.5 mm FWHM ~ 10/t). Using the equation(60) of Jakeman [2(J for 
Gaussian beams, we estimate that the reduction in the enhancement factor can only be in the range of 0.01 to 0.04 
for the different It used. In spite of our best efforts, higher enhancements were not observed. The CBS peak from a 
piece of white paper with a broad peak (FWHM of lOmRad) was observed to have an enhancement of only about 
1.85 without deconvolution. The exact reason for this lower enhancement is not clear, but it could be due to some 
of the laser light being scattered by the non-ideal optical elements, cuvette walls, aggregates and dust which cause 
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an additional component to the diffuse background. This would, however, be sample- independent. Assuming this 
additional contribution to be uniformly distributed, we have attempted to quantify it by studying the enhancement 
factor with respect to absorption in the sample. We conclude from this procedure that the ratio <Jb/a s < 0.1 and that 
the g-factor of milk is greater than 0.55, which is what we expect from our theoretical considerations given earlier. 

The transport properties can also be obtained from the plot of transmission versus the optical thickness p2j . The 
results of the transmission measurements through milk-water suspension is shown in figure-6. At low densities (re- 
gion A), the transmission falls exponentially depending on the extinction coefficient, the slope of the logarithmic 
transmittance being given by a s + a a . In the multiple scattering domain (region B), the transmission falls inversely 
with concentration or optical thickness. At large optical thickness or concentration (region C), the transmission falls 
exponentially depending on the absorption coefficient. The slope of the logarithmic transmittance is <j a in this region. 
These yield a value of the scattering length (l s ) of 130/xm and an absorption coefficient of 1.6cra -1 for milk at the 
largest concentration used. As it can be noticed the absorption coefficient lies just in the region where the CBS 
profiles become insensitive to absorption. Using a value of l t — 470/im obtained from the CBS measurements for 
this concentration and l s — ISOfim from the transmission experiment, a value of g = 0.74 is obtained which appears 
reasonable for such forward scattering media. The inset in figure-6 shows the angle-resolved scattering data, together 
with a theoretical fit to the Henyey-Greenstein function (solid line) with g as the free parameter to be fitted. The fit 
yielded a value of about 0.70 for the anisotropy factor (g) with an error of about 5%. 



V. CONCLUSIONS 



In conclusion, we have investigated the application of Coherent backscattering to measure the light transport pa- 
rameters in tissues using tissue phantoms (milk and dye doped milk), paying particular attention to the limitations 
of the technique. This method yields very good estimates of the transport length but it cannot be used to estimate 
the inelastic length when the absorption of the sample is small (a < 1.0cm -1 ) for the long transport length typical 
of tissues. The maximum measurable absorption coefficient (a ~ 30cm -1 ) with our setup appears to be limited by 
the sensitivity and noise levels of the detector used. It has been pointed out that in principle this technique can also 
be used to measure the anisotropy factor(g) by estimating the single scattering contribution to the angle independent 
diffuse intensity. The single scattering contribution becomes very small for g>0-Q (~ 2% of the total backscattered 

intensity) and the technique becomes insensitive for larger values of g. Thus it does not seem to be a suitable technique 
to estimate the g-factor of biological media which typically are highly anisotropic (g ~ 0.8 to 0.95). Using the value 
of the transport length obtained from CBS measurements and scattering length from transmission measurements, 
a value of 0.74 has been obtained for the anisotropy parameter of milk. This value is confirmed independently by 
angle-resolved scattering measurement where a value of 0.70 is obtained. 



APPENDIX 

For the convenience of the reader, the result of Ref. plj for the angular shape of the CBS peak from a slab of finite 
thickness is reproduced here. 

Let us define [3 = y/(hk/3 + q^), q± = (27r/A)sin0, 5 = (2tt/A)(1 - cos0), rj = (1 + l/cos0)/2Z t and z = 0.71Z t , 
where 9 is the backscattering angle, A is the wavelength of the light, l t is the mean transport length and ij is the 
inelastic length. Now the angular shape of the CBS peak in our notation is given by 



2J?/3 sinh[/3(L + 2z )] W + 5 2 - (i 2 ) 2 + (2/3S) 2 

{ ^-(/3 2 + S 2 - rf ) sinh[/3(L + 2z )} sinhfaL) + 2(ry 2 + S 2 + f3 2 ) cos(SL) 
V 

+ 2(r] 2 + 5 2 - f3 2 ) cosh[/3(L + 2z Q )} cosh(?/L) + 4/3r? sinh(/3L) sinh^L) 

- 2{r) 2 + S 2 + 1 ) cosh(/3L) cosh^L) - 2(?y 2 + S 2 - (3 2 ) cosh^zo) cos{5L) 

- /3?7sinh(2^z )sin( ( 5i) }, (4) 
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where L is the thickness of the slab. It should be noted that the diffuse part of the backscattered light only 
has a weak kinematic dependence on the angle and the ji(9,L) ~ j c (0,L). In the exact backscattering direction 
7l(0,L)=7c(0,i). 
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FIGURE CAPTIONS 



Figure-1 : Plot of the ratio of the backscattering cross-section^) to the total scattering cross-section(cr s ) as a 
function of the anisotropy factor (g) (normalized to unity for isotropic scattering. These quantities were calculated 
using Mie scattering theory for homogeneous spheres. Figure-la is for monodisperse particles and figure-lb shows the 
dependence averaged over a particle size distribution of 30%. 

Figure-2 : Experimental set-up for observing the coherent backscattered light. L1-L3 : lenses; M : mirror; Al,A2 : 
apertures; PI, P2 : polarizers; BS : beam-splitter; BD : beam-dump; S : sample. 

Figure-3 : Measured CBS intensities for milk of two different concentrations. The solid lines arc the calculated fits 
to the corresponding experimental points. The fitted profiles are for : i) l t of (490 ± 20)/im(») and (1100 ± 50)/im(o) 
ii) h of (15000 ± 5000)^m(«) and (13000 ± 5000)^m(o) iii) 7s of (0.3 ± 0.02) for both fits. The inset shows the 
variation of inverse of transport length with the concentration. 

Figure-4 : CBS profiles from methylene blue doped milk for different absorption coefficients. Solid lines are the 
calculated best fits to the corresponding experimental points with i) l t of (670 ± 30)/im(A) and (680 ± 30)/im(»), ii) 
li > 7000/ito(A) and k of (1900 ± 200)/um(») respectively. The inset shows the CBS profile from dye doped milk at 
higher dye concentration (a <~ 10cm -1 and l t = (630 ± 50)/xm). 

Figure-5 : The CBS intensities of the best-fits in figure-4 for an infinite resolution. The enhanced weight of the 
single scattering contribution for the sample with higher absorption resulting in lower enhancement is clearly seen. 

Figure-6 : Variation of the transmittance of milk as a function of the effective optical thickness of milk-water solu- 
tion. The exponential fits in regions A and C are shown. The slopes of the logarithmic transmittances are 8.01 and 
0.16 and yield values of 130/um for the scattering length and 1.6cm -1 for the absorption coefficient respectively, at 
the largest concentration of the milk. Inset shows the angle resolved scattering data. The solid line is a theoretical 
fit to the Hcnyey-Greenstein function giving a value of 0.7 ± 0.04 for g. 
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